arXiv:hep-ph/9808385v2 28 May 1999 


CERN-TH/98-261 
BUTP-98/20 
MPI-PhT/98-58 
DESY 98-101 
|hep-ph/9808385 

August 1998 


Next-to-Leading Order QCD Corrections 
to the Lifetime Difference of Bg Mesons 


M. Beneke, G. Buchalla 
Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

C. Greub 

Institut fiir Theoretische Physik, Universitat Bern 
Sidlerstrasse 5, CH-3012 Berne, Switzerland 

A. Lenz 

Max-Planck-Institut fiir Physik - Werner-Heisenberg-Institut 
Fdhringer Ring 6, D-80805 Munich, Germany 

and 

U. Nierste 

DESY - Theory Group, NotkestraBe 85, D-22607 Hamburg, Germany 


Abstract 

We compute the QCD corrections to the decay rate difference in the Bg-Bg system, 
AP, in the next-to-leading logarithmic approximation using the heavy quark expansion 
approach. Going beyond leading order in QCD is essential to obtain a proper matching 
of the Wilson coefficients to the matrix elements of local operators from lattice gauge 
theory. The lifetime difference is rednced considerably at next-to-leading order. We find 
(AP/P)b^ = (/5^/210MeV)^[0.006i?(mb) -|- 0.150Bs{mi,) — 0.063] in terms of the bag 
parameters B, Bg in the NDR scheme. As a farther application of onr analysis we also 
derive the next-to-leading order result for the mixing-induced CP asymmetry in inclusive 
b —>• uud decays, which measures sin 2a. 


PACS nnmbers: 12.38.Bx, 13.25.Hw, 14.40.Nd 





1. Introduction. The width difference (Ar/r)^^ of the Bg meson CP eigenstates [|^ 
is expected to be about 10 — 20%, among the largest rate differences in the 6-hadron 
sector [0, and might be measured in the near future. A measurement of a sizeable 
(Ar/r)^. would open up the possibility of novel CP violation studies with Bg mesons 
ii- In principle, a measured value for AP Bs could also give some information on the 
mass difference AMb^ 0, if the theoretical prediction for the ratio (AP/AM)^^ can 
be sufficiently well controlled p. Furthermore, as pointed out in [0, if non-standard- 
model sources of CP violation are present in the Bg system, AP b^ can be smaller (but 
not larger) than expected in the standard model. For this reason a lower bound on the 
standard model prediction is of special interest. 

The calculation of inclusive non-leptonic 6-hadron decay observables, such as AT b^ , 
uses the heavy quark expansion (HQE). The decay width difference is expanded in 
powers of AqcB/mb, each term being multiplied by a series of radiative corrections in 
ag{mb). In the case of (AF/F)^^, the leading contribution is parametrically of order 
{Aqcb/' mb)^■ In the framework of the HQE the main ingredients for a reliable pre¬ 
diction (less than 10% uncertainty) are a) subleading corrections in the 1/mb expansion, 
b) the non-perturbative matrix elements of local four-quark operators between H-meson 
states and c) 0{ag) radiative corrections to the Wilson coefficients of these operators. 
The hrst issue has been addressed in P]. The hadronic matrix elements can be studied 
using numerical simulations in lattice QCD. In this letter, we present the next-to-leading 
order QCD radiative corrections to the Wilson coefficient functions for AT b^ ■ In addi¬ 
tion to removing another item from the above list and reducing certain renormalization 
scale ambiguities of the leading order prediction, the inclusion of 0{ag) corrections is 
necessary for a satisfactory matching of the Wilson coefficients to the matrix elements 
to be obtained from lattice calculations. 

Our results provide the hrst calculation of perturbative QCD effects beyond the lead¬ 
ing logarithmic approximation to spectator effects in the HQE for heavy hadron decays. 
The consideration of subleading QCD radiative effects has implications of conceptual 
interest for the construction of the HQE. Soft gluon emission from the spectator s quark 
in the Bg meson leads to power-like IR singularities in individual contributions, which 
would apparently impede the HQE construction, because they cannot be absorbed into 
matrix elements of local operators. It has already been explained in |^, how these severe 
IR divergences cancel in the sum over all cuts of a given diagram, so that the Wilson 
coefficients of four-quark operators relevant to lifetime differences, such as AT Bs , are 
free of infrared singularities. This infrared cancellation is conhrmed by the result of our 
explicit calculation. 

Using the HQE to hnite order in AqaB/'^h rests on the assumption of local quark- 
hadron duality. Little is known in QCD about the actual numerical size of duality- 
violating effects. (See, e.g., Q for a recent discussion of the issue.) Experimentally no 
violation of local quark-hadron duality in inclusive observables of the R-meson sector 
has been established so far. In 13 it has been shown that for AF^^ local duality holds 
exactly in the simultaneous limits of small velocity {Aqcd "C — 2mc -C mb) and large 
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number of colours {Nc —>■ oo). In this case 


Ar\ „ 


l2-A^TB^ 

mb 


0.18. 


( 1 ) 


It is interesting that the numerical value implied by the limiting formula (|^ appears 
to be quite realistic and is in fact consistent with the results of more complete analyses 
1^. The duality assumption can in principle be tested by a confrontation of theo¬ 


retical predictions, based on the HQE, with experiment. This aspect is another major 
motivation for computing AT^,, accurately. 

It is clear from these remarks that a detailed theoretical analysis of AT, and par¬ 
ticularly of 0{as) corrections, is very desirable, both for phenomenological and for con¬ 
ceptual reasons. In this letter we shall concentrate on the presentation of our results 
and a brief discussion of their main aspects. Details and an extension of our analysis to 
other fe-hadron lifetime differences will be given in a forthcoming publication 0- 


2. Formalism and next-to-leading order results. In the limit of CP conservation 
the mass eigenstates of the Bg-Bg system are \Bh/l) = (l-Bs) ± l-Ss))/'\/2, using the 
convention CP\Bg) = —\Bs). The width difference between mass eigenstates is then 
given by 

ATs^ = Ti — T//= —2 ri2 = —2 r2i, (2) 

where Tij are the elements of the decay-width matrix, i,j = 1, 2 (|1) = \Bg), |2) = \Bg)). 
In writing we assumed standard CKM phase conventions |T^. For more information 
about the basic formulas see for instance |p (and references therein). 

The decay width is related to the absorptive part of the forward scattering amplitude 
via the optical theorem |Q. The off-diagonal element of the decay-width matrix may thus 
be written as 

r2i =- {Bg\T\Bg). (3) 

2Mb/ "II"/ V ; 

The normalization of states is {Bg\Bg) = 2EV (conventional relativistic normalization) 
and the transition operator T is dehned by P] 

T = (4) 

Here Tfg// is the low energy effective weak Hamiltonian mediating bottom quark decay. 
The component that is relevant for r 2 i reads explicitly 



W.// = 1 

^^j^CrQr + CsQsy 

(5) 

with the operators 




Qi = 

(bi Cj ) Y ( Cj Si^y 


(6) 

Q 3 = 

{biSi)v-A{QjQj)v-A 

Qa {pi^j^V—A{,QjQi)v—A’! 

(7) 
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Qs —A{^Q.jQ.j)v+A Qq j^V—A{^Q.jQ.i^V-\-A^ 


( 8 ) 

(9) 


Here the z, j are colour indices and a summation over q = u, d, s, c, b is implied. V i: A 
refers to 7'^(1 ± 75 ) and S — P (which we need below) to (1 — 75 ). Ci,... ,Cq are the 
corresponding Wilson coefficient functions, which are known at next-to-leading order. 
We have also included the chromomagnetic operator Qs, contributing to T at 0{as). 
(Note that for a negative Cg, as conventionally used in the literature, the Feynman rule 
for the quark-gluon vertex is A detailed review and explicit expressions may 

be found in [^. Cabibbo suppressed channels have been neglected in (|^). 

Expanding the operator product (Q) for small x ~ l/m;,, the transition operator T 
can be written, to leading order in the l/rrih expansion, as 


T = 


Glml 

127r 


(V;K.)" [F(z)Q(ti2) + Fs(z)Qs{p,)] 


with 2 ; = ml/ml and the basis of AB = 2 operators 

Q = {biSi)v-A{bjSj)v-A, Qs = (biSi)s-p{bjSj)s-i 


( 10 ) 


( 11 ) 


In writing (^) we have used Fierz identities and the equations of motion to eliminate 
the colour re-arranged operators 


Q {piS j^v—A{bjS/)v—At Qs j) S—piP j^i) S—I 


( 12 ) 


always working to leading order in 1/mb (see below). The Wilson coefficients F and 
Fs can be extracted by computing the matrix elements between quark states of T in 
(I) (Fig- 0), as well as those of Q and Qs, to a given order in QCD perturbation 
theory, and comparing them with (W). (The external b quarks are taken to be on- 
shell.) This matching procedure factorizes the perturbatively calculable short-distance 
contributions (Wilson coefficients) from the long-distance dynamics, parametrized by the 
non-perturbative matrix elements of local AB = 2 operators. We do not use heavy quark 
effective theory (HQET) to expand these matrix elements explicitly in l/mt- They are 
to be understood in full QCD, in accordance with the treatment of 1/mb effects in [^. 

In our discussion we shall first concentrate on the contribution to T from the operators 
Qi and Q 2 in (|^). The penguin operators have small coefficients and are numerically 
less important. Their effect will be included later on. Thus for the time being only the 
diagrams Di - Dio in Fig. Q are considered, while Du and D 12 belong to the penguin 
sector to be discussed below. 

Working in leading order, the matching calculation for ([T0| ) has to be performed to 
zeroth order in as- At next-to-leading order the coefficients C' 1,2 have to be taken in next- 
to-leading logarithmic approximation |]T^ |T^ and the 0{as) matching corrections to the 
coefficients in m have to be computed. At this order the coefficients in (pUj ) depend 
on the renormalization scheme chosen for their evaluation. This scheme dependence is 
cancelled by the matrix elements of the operators in ([T(]|), which have to be determined 
accordingly. 
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Figure 1: 0{as) corrections to the Bg-Bg transition operator, Di - Di 2 - Also shown are 
the corrections to the matrix elements of local AB = 2 operators, Ei - i? 4 , required for 
a proper factorization of short-distance and long-distance contributions. Not displayed 
explicitly are E 2 , D[, D^i D^, Dj, Dg, D[q and 0 ^2, which are obtained by 
rotating the corresponding diagrams by 180°. 


The scheme we employ for our result is specihed as follows. We use dimensional regu¬ 
larization with anti-commuting 75 and modihed minimal {MS) subtraction of ultraviolet 
singularities (NDR scheme). In addition we project D-dimensional Dirac-structures, aris¬ 
ing at intermediate stages of the calculation, according to the prescriptions (H = 4 — 2e) 

[7^7“7‘^(1 - - l5)]ki ^ (16 - M) [ 7^(1 - 75)]7 [ 7^(1 - 75 )]^^, (13) 


- 75)]p [ 7 ;^ 7 « 7 /.(l - 75)]fc/ ^ (4 - 8£)[7^(1 - 75)]p- [7^,(1 - 75)] 


kl, 


(14) 


[ 7 /^ 7 ^( 1 - 75 )].^. ^ (8-4£)[l-75]i^- [i-^5]^;-(8-8£)[l-75]iz [l-75]fci- (15) 

The projections are equivalent to the subtraction of evanescent operators, dehned by 


the difference of the left- and right-hand sides of (|T3|-^). The dehnition of evanescent 
operators is discussed in great detail in [|^ |^, ^ and we use the basis of 


171 for 


the projection. The prescriptions (p!3|-[T5|) complete the dehnition of our renormalization 
scheme. They preserve Fierz symmetry, i.e. the one-loop matrix elements of Q, Qs and 
Qs equal the matrix elements of the operators that one obtains from Q, Qs, Qs 
(4-dimensional) Fierz transformations. The projections (^-K) are sufficient if we use 
the Fierz form of Qi, Q 2 in (^ that corresponds to closed charm-quark loops in Fig. |I[ 
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Since the renormalization of Qi, Q 2 respects Fierz symmetry, this choice can always be 
made. 


We now give the result for the transition operator (W) at next-to-leading order, still 
neglecting the penguin sector. The coefficients in m can be written as 


F{z) — + T\ 2 (^)C'l(/^l)C' 2 (hl) + -P 22 (^)C '2 (hi)) 


F,,(z) = F}f>iz) + 




dvr 


and similarly for Fs{z). The leading order functions read explicitly 

fI°\z) = (1 - z) F^%{z) = 3v/r^4i (1 + 2z), 

fI^{z) = 2VT^ (1 - F^%(z') = 2VT^ (1 + 2z), 


F^°\z) = ^(l-4zr/^ 


^( 0 ) 

S,22 
( 1 ) l?(l) 


{z) = —\/l — 4z (1 + 2z). 


The next-to-leading order expressions F^j , FgJj are 


fI^P(z) = 32(1 — z)(l — 2z) ^Li 2 ((T^) -|- In^ cr -|- - lncrln(l — 4z) — In a In 2 ;^ 

-h 64(1 - z)(l - 2z) ^Li 2 (cr) -h ^ ln(l - a) Ina^ 

- 4(13 - 26z - 4^2 -1- 14z^) In a 


+ Vl -4^ 4(13 - 10^) Inz- 12(3 - 2z) ln(l - 4z) 


+-(109 - 226z + 168^2)1 + 2Vl-4z(5 - 8z) In 
6 J mb 


4!n(^) 


= 32(1 — 4z‘^) ^Li 2 ((T^) + In^ ^ ln(Tln(l — 4z) — In cr In z 

+ 64(1 — 4z‘^) ^Li 2 ((T) + - ln(l — a) In cr^ 

- 16{4-2z-7z^+ 14z^)lna 
+ Vl — 4;^ 64(1 + 2z) Inz — 48(1 + 2z) ln(l — 4z) 


--(1 - 6z)(5 + 7^)1 - 32V1 -4^(1 + 2z) In 
3 J mb 


-^ 12 ^ (^) = “^(1 ~ ^)(1 ~ 2z) ^Li 2 (cr^) + In^ cr + - ln(Tln(l — 4z) — In cr In 

+ ^(1 - 2)(1 - 2s) +2(17) + ^ln(l - CT)liih 


(16) 

(17) 

(18) 

(19) 

( 20 ) 


( 21 ) 


( 22 ) 


5 


















+ (2 - 259;^ + 662^2 - 7Qz^ - 200z'^) 


Inz 


In a 

6z 


— \/i — 4z 


(2 - 255z + 316z^) -+ 8(3 - 2z) ln(l - Az) 


6z 


+ -{127 - 199z - 75z^) 
9 


2VT^^{17 - 26z) In — + -Vl - 4z(5 - 8z) In —, 

mb 3 mb 


+2 


(23) 


= “^(1 ~ 42:^) fLi 2 ((J^) + In^ a + - ln(Tln(l — Az) — In a In 2 : 
o \ Z 

+ ^(1-4^^) (^Li 2 (cT) + ^ln(l - cT)lna^ 


+ (1 - 35;^ + 4^2 + 76z^ - lOO/) 


4 Ina 
3z 


— \/l — 4z 
4 


(1 - 33z - + 32(1 + 2z) ln(l - 4z) 


3z 
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(68 + 49z - 150z" 


I 6 VI — 4^(1 + 2z) In —— — Vl — 42;(1 + 2z) In —, 

mb 3 mb 


(24) 


4‘’(z) = 


(4 — 2I2: + 2 z^) {hi2{(J^) + In^ a + - In a ln(l — Az) — In a In 2: 
+ ^(1 - 22:) (5 - 2 V ^Li 2 (a) + ^ ln(l - a) Ina^ 


(7 + 13z - 194V + 304V - 64V)— - —(1 - 10V 

Qz 3 

(7 + 27z - 250V)^ - 4(1 - 6z) ln(l - Az) 


+ Vl — 4;^ 


6z 


— (115 + 632^ + 96V) 
18 

+1 


— 2 VI — 42:(5 — 2z) In — H—Vl — 42:(2 — 5z) In —, 

mb 3 mb 

32 / I \ 

Fg^22{^) = f Li2(V) + iV a + - ln(Tln(l — 4 z) — Xnalnzj 

+ y (1 - 4V) (^Li 2 (a) + ^ ln(l - a) \na^ 


(25) 
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+ (1 + 7^ + lOz^ - Q8z^ + + —(1 + 2z) 

3z 3 


a/I — Az 


A l-p 2^ 

(1 + 92; + 262 ;^)-16(1 + 22 ;) ln(l — Az) 


3z 


+-(19 + 53^ + 24^2) 

y 

— 16 a /1 — 42;(1 + 22 ;) In — + — Vl — 4^(1 + 22 ;) In —. 

mf, 3 rUb 


In these equations we have set Nc = 3 and used 

t) 


(26) 


Li 2 (a;) = 


/o t 


a = 


1 - ^/l-Az 
1 + \/l — Az 


m: 


2: = 


mt 


(27) 


The dependence on the renormalization scale pi in (pl| -p6D cancels against the scale 
dependence of the Ai? = 1 Wilson coefficients C'/pi) in (PT|) to the considered order in 
ttg. Likewise, the dependence on /i 2 is cancelled by the matrix elements of the Ai? = 2 
operators Q and Qs- To check this, we note that the scale dependence of matrix elements 
(Q) = ((Q), {Qs), {Qs)Y is given by 


where 


= -/7f «), 

d In 47r 

(28) 

0 

0 


7^°) = 0 -28/3 4/3 . 

(29) 


Vo 16/3 32/3, 


The operator Qs, which is redundant at leading power in l/rrib, can then be eliminated 

/i 2 -independence can be verified. 

correspond to the use of the one-loop pole mass for the b quark 


using (|35D below and the /i 2 -independence can be verified. 
Our results in (p!- 


in (p!0|) . In 2 ; = 777 . 2 / 777.2 there is no difference between the ratios of pole masses and MS 
masses at next-to-leading order in ag. 

We next present the contribution from QCD penguins to the transition operator T. 
This sector has been treated in 0 in the leading logarithmic approximation. At next- 
to-leading order the penguin coefficients C^,, ■ ■ ■ ,Cq have to be computed with next-to- 
leading logarithmic accuracy JT^, |T^ and the diagram Du in Fig. |l| must be evaluated. 
Because the coefficients O 3 ,... ,Cq are small (~ few per cent), contributions of second 
order in these coefficients are safely negligible and it is sufficient to calculate only the 
interference of Cs,..., Ce with Ci and +' 2 - A consistent way to implement this approxi¬ 
mation at NLO is to treat C ^,..., Cg as formally of 0{as). The standard NLO formula 


(as reviewed in |T^) can be used for the penguin coefficients, except that terms of order 
OgCs ,..., ttsCg have to be dropped. Accordingly, only current-current operators Qi, Q 2 
are inserted into the diagrams Di - Du. The NLO result thus obtained is manifestly 
scheme independent and formally of order 0 ( 03 ^...^g) = 0{as). A further contribution. 
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absent in leading logarithmic approximation, comes from the chromomagnetic operator 
Qs in @ and is shown as diagram D 12 of Fig. |l|. Since this contribution arises hrst at 
0{as), the lowest order expression is sufficient for Cs- For the NLO penguin contribution 
we then hnd 

T, = . (30) 

P(z) = VT^ ((1 - + 1(1 - 4z)R-M + SzRKjz,)'^ + ^4^F,(z)CI(mi), 

(31) 

Ps(z) = ^r^(l + 2z) (K'M - KM) - (32) 


Fp{z) = --Vl - 42:(1 + 2z) 


(33) 


/ii 


3C'8(/ii) 


21n — + - + 4^ - ln ;2 + a/I - 42:(1 + 2z) In cr + , 

mb 3 C2(/ii) 


where we dehned the combinations K[ = 2{?)CiC^ + C'iC '4 + C' 2 C' 3 ), = 2 C 2 C 4 and 

/F 3 = 2 ( 3 C'iC' 5 +CiC' 6 +C 2 C 5 +C' 2 C 6 ). The explicit /xi-dependence in Fp[z) is cancelled by 
the /ii-dependence of (Fs, ..., Cq. Note that for the penguin sector the scale and scheme 
dependence of the matrix elements of Q and Qs are effects beyond the considered order 
and numerically negligible. 

Beyond leading order in the 1/mb expansion, several further operators contribute to 
AF, denoted by Ri in |p. Inspection of the factorized matrix elements of these oper¬ 
ators given in shows that these superhcially power-suppressed operators contribute 
at leading power beyond tree level. This is due to the fact that the operators are de¬ 
hned in QCD (rather than in HQET) and hence the one-loop matrix elements contain 
mb explicitly. The leading power piece arises only from loop momenta of order mb and 
can therefore be subtracted perturbatively. This subtraction is necessary for a complete 
calculation of the as correction at leading power and is taken into account in the above 
result. To be more explicit, consider as an example the operator 

Ro = Qs + Qs + -^Qi (34) 

which can be reduced to an explicitly power-suppressed operator using Fierz transfor¬ 
mations and the equations of motion. (We have used this in our calculation to eliminate 
Qs-) At order Os, we hnd the leading power contribution 


iK) = 


dvr 




2 -Nr 


(W+l)ln^ + 

mb N, 

h 


(Q) 


(35) 


+ 4(Ac -|- 1 ) In-h 2{Nc + 1 )) {Qs) 

mb 


+ 0{l/mb). 


It is crucial that this relation holds independent of the external state, so that power 
counting is again manifest after subtracting the right hand side of (^^ from the matrix 
















element of Rq. (The procedure discussed here bears some similarity with mixing of higher 
dimension operators into lower dimension operators in cut-off or lattice regularizations.) 
These subtractions must be kept in mind when a non-perturbative evaluation of the 
matrix elements of the Ri is combined with the present NLO results. In the factorization 
approximation of 0 these subtractions correspond to using the pole b quark mass in the 
expressions for the factorized matrix elements. Indeed, in the > oo limit, we hnd 0 


{Bs\Rq\Bs) 




Ml 


{mb + m, 


^ \2 


- f. 


I Ml 

IDs 


dsNc 

47r 



(36) 


up to corrections of order {Ncttg 
expression of 


and hoicD/mb. The 6-quark mass rfit in the second 
is the MS mass at the scale fi, which corresponds to our renormal¬ 
ization of the scalar operators Qs, Qs- To obtain the third expression we used the 
1 -loop relation between the pole and MS mass in the large-A^c limit and the fact that 
Mbs ~ ^b,poie = C>(Aqcd)- The same result as ( |36[) is obtained from the large-limit 
of (1351) . (In deriving (|36|) we have used the Fierz transform of Qs- For the coincidence 
of (^) with (|35|) it is crucial that the choice in ([T5|) maintains the Fierz symmetry in 
the one-loop matrix elements entering 


3. Discussion. The complete expression for AT Bs with short-distance coefficients at 
NLO in QCD is given by 

AF^ 167i^B{Bs ^ Xeu) ^Mb^ 

r )b. 9{z)vqcd ml ' ^ ^ 

where 

G{z) = F{z) + P{z) and Gs{z) =-{Fs{z) + Ps{z)). (38) 

We eliminated the total decay rate F b^ in favour of the semileptonic branching ratio 
B{Bs —>■ Xeu), as done in 0. This cancels the dependence of (AF/F) on Wb and 
introduces the phase space function 

g{z) = l-8z + 8z^-z^-12z‘^\nz, (39) 


as well as the QCD correction factor 

2as{mb) 


Vqcd = 1 


37r 




TT 


V-z? + I 


The latter is written here in the approximate form of ^1 
parametrize the matrix elements of Q and Qs, 


(40) 

The bag factors B and Bs 


{Bs\Q\Bs) 

{Bs\Qs\Bs) 


flMlB, 


- If ^ 




9 Ml 

Ml TZ - 

“ {mb -I- m^)^ 


Be 


(41) 

(42) 
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hi 

mb/2 

mb 

2mb 

Gs = -{Fs + Ps) 

0.743 

0.937 

1.018 

G!?' = -(Ff + Pf) 

1.622 

1.440 

1.292 

G = F + P 

0.023 

0.030 

0.036 

q { 0 ) ^ ^(0) _|_ p(0) 

0.013 

0.047 

0.097 

-Fs 

0.867 

1.045 

1.111 

p(o) 

-^S 

1.729 

1.513 

1.341 

F 

0.042 

0.045 

0.049 

_p(0) 

0.030 

0.057 

0.103 


Table 1: Numerical values of the Wilson coefficients G, Gs, F, Fs, at next-to-leading 
order in the NDR-scheme with evanescent operators subtracted as described in the text 
(for fi 2 = mb). Leading order results (superscript ‘(0)’) are also shown for comparison. 


The masses mb = mb{fi), ms refer to the MS dehnition. The relation of mb to the pole 
mass mb is 


mfc(/i) = rUb 



(43) 


Finally, 6i/m describes l/mb corrections. The explicit expression for 6i/m in the factoriza¬ 
tion approximation can be found in |^. In the present NLO approximation the 6-quark 
mass that appears in 6i/m is the pole mass as mentioned above. 

For the numerical evaluation we use the following input parameters (central values): 


mb = 4.8 GeV mb{mb) = 4.4 GeV 2 ; = 0.085 mg = 0.2 GeV, (44) 

Mb, = 5.37 GeV fs, = 0.21 GeV B{Bs Xeu) = 0.104. (45) 

The two-loop expression is used throughout for the QGD coupling ag in the form given 
in []^ with = 0.225 GeV. The NLO coefficients F{z), Fs{z) in ( |TI]|) and P{z), 
Ps{z) in (^) are consistently expanded to hrst order in ag. We take fii = = mb as 

central values for the renormalization scales. The dependence of AF b, on rhg is marginal 
and its dependence on 2 ; stems almost totally from g(z) in (^). 

The results for the Wilson coefficients are displayed in Table |^. The contribution 
of Qs dominates AF b, since the coefficient of Qs is numerically much larger than the 
one of Q. These coefficients are independent of the scale fii, related to the AB = 1 
operators, up to terms of next-to-next-to-leading order. As can be seen from Table ^ 
the residual /ii dependence is indeed substantially decreased for the coefficient G of Q. 
The reduction of scale dependence is less pronounced for the coefficient Gs of Qs- On 
the other hand, in our renormalization scheme, the central value at NLO is considerably 
smaller (by about 30%) than the leading order result. However, due to the scheme 
dependence of the Wilson coefficient, it is premature to draw dehnitive conclusions on 
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the size of AFwithout combining the coefficient functions with B and Bs, computed 
in the same scheme. For quick reference we rewrite (|^) as 


AFX 

~)b.. 




210 MeV 


[0.006 5(mfe) + 0.150 Bs{mh) - 0.063], 


(46) 


where the numbers are obtained with our central parameter set. 

A preliminary lattice evaluation of the relevant bag parameters can be found in 
together with a complete 1-loop lattice-to-continuum matching. We can use their result 
to obtain a (conceptually) complete (but numerically preliminary) next-to-leading-order 
result for (AF/F)^^. In a continuum renormalization scheme different from ours for 
the operators Qs and Qs is chosen. We computed the relation between the two schemes 
to 0{as) and found 


Bs 

Bs 


Bt 


Ots 

dvr 

(y.Q 




B+ + + -B-* 

^ 47r I 6 * ^ 2 " 


(47) 

(48) 


Here Bs is the bag parameter related to Qs, dehned in analogy with (|4^) with Qs —>■ 
Qs, Bs —> Bs, and the numerical factor (—5/3) —(-1-1/3). Bf and B^ are the bag 
parameters Bs and Bs, respectively, but in the scheme of |^. From the quoted estimates 
Bt ifJ-o) = 0.80 ± 0.01, H^(/io) = 0.94 ± 0.01 (where the errors are statistical only 
and /io = 2.33 GeV) we infer BsifJ^o) = 0.81, BsifJ^o) = 0.87. Using (P^), and taking the 
running of fhh (^) into account, the bag parameter at the scale mb is given by 

Bsimb) = Bsifio) + ^ In ^ = 0.75. (49) 

dvr po V 3 15 / 

We also take B{mb) = 0.9 from the compilation ||^. Without the 1/mb corrections, 
(1^) then becomes (AF/F)^^ = 0.118 (/b^/ 210 MeV)^. There is a ±15% error from 
(continuum) perturbation theory, obtained from varying pi between mb/2 and 2mb, and 
a negligible ±1% statistical error from the lattice simulation. In addition the sizeable 
1/mb correction, computed in 0, has to be included. The estimate —0.063 in (|46D is 
obtained using factorization of hadronic matrix elements and has a relative error of at 
least ±20%. As a preliminary result we may therefore write 

= ( 210^Mev ) (0-054lo;o32 (/^i-dep.)±??? (latt. syst.)) . (50) 


We emphasize the preliminary nature of the central value which depends crucially on 
the estimate Bt{^o) = 0.80 taken from There is an unspecified systematic error 


(indicated by the question marks) attached to this number, related to the fact that the 
lattice calculation has been performed in the quenched approximation, at hnite lattice 
spacing and with a “6-quark” mass in the charm quark mass region, without extrapo¬ 
lation to the continuum limit and to realistic 6 quark masses, respectively. Clearly, for 
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further progress improved lattice determinations of bag parameters, most importantly 
of Bg, are mandatory, and the numerical value for (Ah/r)^^ above has to be regarded 
in this context. 

The rather low number for the central value, compared to the leading order analysis 
of 0, is a consequence of the fact that l/rnt effects and penguin contributions 0, as 
well as NLO QCD corrections, all lead to a reduction of AT Bs- This is further reinforced 
by the small value = 0.75 in our example. 

It is interesting to consider the ratio of AT to the mass difference AMb^ @, |^, 
in which the dependence on the decay constant cancels out and the sensitivity to 
Vcb is considerably reduced. In addition (AT/AM) Bs only depends on the ratio of bag 
parameters. Generalizing the results given in 0 to include the next-to-leading order 
QCD corrections we can write 


/ AP ^ 

TT ml 

PcfePcS 

2 1 

VAMy 

’ b . 2 Ml 

VtsVtb 

VBSoiXt) 


(51) 


G{z) -F 


Ml 


3 (m^ -|- rris 


G. 


^ -F Vl - 42; 6i/rr}j , 


where t]b is the (scheme-dependent) next-to-leading order QCD factor entering AMb^ 
In the usual NDR scheme riB{mb) = 0.846. The top-quark mass dependent function 
So{xt) = So{{fht/Mw)l = 2.41 for rht = 167 GeV. In analogy with (|50|) we then have 


AT \ 

am) 


B, 


^2.63](°;3g (/ii-dep.) ± ??? (latt. syst.)^ 
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-3 


(52) 


The next-to-leading order calculation presented in this article can also be applied, by 
taking the limit 2; —> 0, to the mixing-induced CP asymmetry in inclusive Bd{Bd) —>■ uudd 
decays The time-dependent asymmetry is given by A{t) = Im,^sin AMt, where the 


coefficient Im^ is a measure of the CKM parameter sin 2a. With next-to-leading order 
QCD corrections included, the expression for Im,^ in Eq. (25) of is modihed to read 


Im^ = —0.12 sin 2a I 


fs 


2 r 


180 MeV 


O.UB + OMBs - 0.07 - 0.06 


sinasin(a -|- /?) 
sin fd sin 2a 


(53) 

In this equation the bag factors B and Bg (taken at p = mb) are the Bd analogues of 
those dehned in (^Tj) and (^2[) for Bg. A detailed discussion of CP asymmetries will be 
presented in 
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